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Abstract 

The phase structure and equation of state for two-flavor quark matter under 
compact star constraints is studied within a nonlocal chiral quark model. Chiral 
symmetry breaking leads to rather large, density dependent quark masses at 
the phase transition to quark matter. The influence of diquark pairing gaps and 
quark masses on density dependent emissivities for the direct URCA is discussed. 
Since m„ > to^, the direct URCA process due to quark masses cannot occur. 
We present cooling curves for model quark stars and discuss their relation to 
observational data. 



1 Introduction 

Emissivities and mean free paths of photons and neutrinos are essential for 
most astrophysical phenomena (supernovae, neutron star cooling, gamma ray 
bursts(GRB), pulsar kicks [5] etc.). Their correct treatment is one of the chal- 
lenging tasks in astrophysics. Theoretical predictions of quark matter prop- 
erties inside compact stellar objects including the possibility of different color- 
superconducting phases, resulted in a number of emissivity calculations in quark 
matter until now (e.g. [D El [3l H] ) . The first calculation of neutrino emissivi- 
ties in quark matter has been done by Iwamoto [T]. He found that the matrix 
element of the direct URCA process would vanish and this important cooling 
process in quark stars could not occur if one neglects quark-quark interactions 
and quark masses. Most calculations consider the effect of quark-quark inter- 
actions to obtain a finite matrix element since the current up and down quark 
masses are small and their influence on the emissivity is negligible. However, 
the up and down quark masses can be up to almost two orders of magnitude 
larger than the current quark mass in the vicinity of the phase transition to 
quark matter. This can be shown within NJL-type chiral quark models [1]. 
The influence of both quark masses and diquark pairing on the direct URCA 
emissivities is discussed in the following sections. 



2 Relativist ic chiral quark model 

The grand canonical potential for quark matter in a 2SC superconductor is 
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with f^e ^ -^yi27r2 - (J-qT'^/G - 7t:^T^/180 being the thermodynamic po- 
tential of ultra-relativistic electrons, where fiQ — —fie, and f2o is the divergent 



vacuum contribution which has to be subtracted to assure vanishing energy and 
pressure of the vacuum. Out of the twelve eigenvalues Aq four belong to the 
ungapped blue quarks and can be determined easily by using textbook methods 
[6] as A1..4 = E-f{p) ± fifb, with the dispersion relation Ef{p) = ^p'^ + M'j (p) 

containing the dynamical mass function Mf{p) ~ mf+g{p)(j)f for the two flavors 
f — u,d. We have introduced the chemical potentials for the quarks of unpaired 
color ^ub = Mb/3 + 2^q/3 - 2^8/2 and fidb = fJ-ut - fJ-Q- The other eight eigen- 
values belong to the red and green quarks which are paired in the 2SC state 
and have therfore identical eigenvalue spectra. It is thus sufficient to determine 
the four eigenvalues for the red quarks by solving a quartic equation similar to 
the one discussed in [7] for the CFL phase. The parameter a describes flavor 
mixing due to instanton induced interactions, see [9l [5] . For the formfactor we 
used a simple NJL cut off g{p) = 9{1 — p/A). The numerical solutions of this 
section are obtained with the parameters from Tabic 5.2 of 5 ; current quark 
mass mu,d = 5.5 MeV , coupHng constant G^A^ = 2.319, cut oflf A = 602.3 MeV 
and diquark coupling Go — VdGs with rjjj — I (strong coupling). Fig.[T]shows 
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Figure 1: Quark masses for asymmetric matter under charge conservation without 
diquark gap and flavor mixing (left panel), with diquark gap and without flavor mixing 
(center), and with diquark gap and flavor mixing (right panel). 

that after the chiral phase transition is first order, signalled by a sudden drop 
of the constituent quark masses at the critical chemical potential /i"'* — 340 
MeV. The region just above the transition, /x^"^'* < /i < 400 MeV, relevant for 
quark stars (or quark cores of neutron stars), shows still nonperturbatively large 
quark masses strongly decreasing towards their current values with increasing 
fi. In asymmetric quark matter under charge conservation the up quark mass is 
always larger than the down quark mass. The magnitude of this mass difference 
decreases with the flavor mixing factor a, being maximal for a — 0. Note that 
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the diquark gap A has a similar but smaller mixing effect on the quark flavors. 



3 Direct URCA emissivities 

The resulting neutrino emissivity from the direct URCA reactions 

d u + + , u + ^ d + i^e 

has been calculated in [T] and can be expressed as 
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Here 6c is the Cabibbo angle, G the weak coupling constant and O^e is the 
angle between the up-quark and electron momenta obtained from the con- 
dition of momentum conservation in the matrix element, see Fig. [2] . We 
take into account that in the 2SC phase 2/3 of the quarks are paired and 
thus the corresponding direct URCA emissivity is suppressed by a factor rj — 
(C(C) + exp[(A - 0.5/ie)/T])"^ taken from Ref. 0, where we choose C(C) = 1. 
For the late cooling stage with temperatures T <C 1 MeV, neutrinos are un- 
trapped and their effects on momentum conservation and chemical potentials 
can be neglected. Trigonometric relations can be used to find an appropriate 




Figure 2: Triangle of momentum conservation for the URCA process 



analytical expression for the momentum conservation from Fig. [2 which in 
lowest order of Ode is given by 
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For small angles holds 9de — due, so that an expression for the matrix element 
of the direct URCA process can be derived. Following Iwamoto [T] we take into 
account either quark-quark interactions to lowest order in the strong coupling 
constant as (O or the effect of finite masses ([6]) 
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From Eqs. (H])-® with the /3-equihbrium condition = ^-u^ Me an expression 
for the angle Ode determining the emissivity Eq. ^ can be found similar terms 
for Eq. Q 




The mass effect for current quark masses is negligible in comparison to the 




Figure 3: Emissivities calculated for parameters obtained from an ideal Fermi gas 
model with the quark current masses rau = 5 MeV and rud = 9 MeV (left panel), 
vs. a relativistic quark model under /3-equilibrium and charge conservation without 
flavor mixing (center) where rriu = 98 MeV and rrid = 45 MeV and with flavor mixing 
a = 0.11 (right panel), where m„ = 62 MeV and md ~ 46 MeV. The diquark Gap 
is in all flgures A ~ 192 MeV. The solid curves correspond to emissivities via quark- 
quark interaction, dashed curves to the mass effect and dash-dotted curves show the 
influence of a diquark gap. The baryon chemical potential is fixed to 1.1 GeV. 

perturbative quark-quark interaction effect, as stated in [1] (see left panel of 
Fig. [3|). However, this statement is restricted to the region of the QCD phase 
diagram, where the perturbative treatment is possible at all and the quark 
masses are of the order of the current values rriu.d ~ 5 — 9 MeV. If > md the 
momentum conservation is violated and the URCA process via the mass effect 
does not work, despite the large quark masses. The influence of the quark masses 
on the Fermi momentum of up and down quarks in a relativistic chiral quark 
model together with the perturbative one-gluon exchange (OGE) interactions 
results in slightly higher emissivities compared to the purely perturbative ones. 
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Emissivities in superconducting quark matter are suppressed by the factor 77, 
leading to lower emissivities, cf. Fig. [31 

4 Cooling 

In order to qualitatively estimate the consequences of the emissivity effects dis- 
cusses above on the cooling behavior of compact stars, we consider the simplified 
model of a homogeneous quark star. The temperature-age behavior is obtained 
by inverting the solution of the cooling equation 

Tit) 

X] 'l^Jt) * " quark, e", 7, gluon; j = v, 7. (8) 

T(to) 

The contributions of the different species to the specific heat and the photon 
contribution to the star luminosity can be found in [TU] . The double logarithmic 




log,o(t[yr]) 



Figure 4: Cooling curves for the corresponding emissivities from Fig. |3]without flavor 
mixing (center) for a homogeneous quark star {/ib = 1.1 GeV) with a radius of 9 km 
and a mass 1.2 Mq. . 

plot Fig. [4] shows the surface temperature Tg of the star vs. its age. The 
relation between crust and surface temperatures is taken from the Tsuruta law 
Ts = (10 Tc)^/^. The cooling curves are sensitive to the quark masses and 
the diquark gap supression. Note that the configuration of a homogeneous 
quark star without hadronic shell is quite unrealistic and the comparison with 
observational data is only made to give an orientation about the relations in 
a T-t plot. A realisitic treatment of quark star cooling including heat flux, as 
well as neutrino and photon transport with temperature and density profiles 
can alter these results strongly. 
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5 Conclusion 



Within our study of the classical URCA process due to perturbative OGE inter- 
actions wc find that constituent quark masses and chemical potentials obtained 
from a NJL type relativistic chiral quark model lead to emissivities comparable 
to the results from the ideal Fermi gas. Neutrino emissivities due to the quark 
masses can only occur if rn,i > m„. Since in the chiral quark model under charge 
conservation the selfconsistent masses obey m<j < m„, momentum conservation 
is violated and the URCA process can not occur. In hybrid stars (neutron stars 
with a quark matter core) or model quark stars we can proceed from the assump- 
tion that the quark matter close to the phase transition is in a non-perturbative 
density region rather than in a perturbative one. It is questionable whether 
OGE interactions are applicable in the non-perturbative density regime close 
to the chiral phase transition. Therefore, non-perturbative effects like, e.g., p 
and u) meson exchanges could be possible candidates to replace the perturbative 
OGE interaction in the non-perturbative density regime. 
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